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The background of gravitational waves produced by the ensemble of rotating neutron stars (which 
includes pulsars, magnetars, and gravitars) is investigated. A formula for fi(/) (a function that is 
commonly used to quantify the background, and is directly related to its energy density) is derived, 
without making the usual assumption that each radiating system evolves on a short time scale com- 
pared to the Hubble time; the time evolution of the systems since their formation until the present 
day is properly taken into account. Moreover, the formula allows one to distinguish the different 
parts of the background: the unresolvable (which forms a stochastic background or confusion noise, 
since the waveforms composing it cannot be either individually observed or subtracted out of the 
data of a detector) and the resolvable. Several estimations of the background are obtained, for 
different assumptions on the parameters that characterize neutron stars and their population. In 
particular, different initial spin period distributions lead to very different results. For one of the 
models, with slow initial spins, the detection of the background by present or planned detectors 
can be rejected. However, other models do predict the detection of the background, that would be 
unresolvable, by the future ground-based gravitational wave detector ET. A robust upper limit for 
the background of rotating neutron stars is obtained; it does not exceed the detection threshold of 
two cross-correlated Advanced LIGO interferometers. If gravitars exist and constitute more than 
a few percent of the neutron star population, then they produce an unresolvable background that 
could be detected by ET. Under the most reasonable assumptions on the parameters characterizing 
a neutron star, the background is too faint to be detected. Previous papers have suggested neutron 
star models in which large magnetic fields (like the ones that characterize magnetars) induce big 
deformations in the star, which produce a stronger emission of gravitational radiation. Considering 
the most optimistic (in terms of the detection of gravitational waves) of these models, an upper 
limit for the background produced by magnetars is obtained; it could be detected by ET, but not by 
BBO or DECIGO. Simple approximate formulas to characterize both the total and the unresolvable 
backgrounds are given for the ensemble of rotating neutron stars, and, for completion, also for the 
ensemble of binary star systems. 



I. INTRODUCTION 

The topic of this paper is the gravitational wave back- 
ground p] [2] produced by the ensemble of rotating neu- 
tron stars in the universe. These systems are modeled 
as isolated neutron stars [3] that are formed with an ini- 
tial spin frequency, and lose energy via electromagnetic 
dipole emission [H [5] and via quadrupolar gravitational 
radiation [SJ |7| . The ensemble of rotating neutron stars 
contains the populations of pulsars, magnetars, and grav- 
itars. 

Pulsars [5] are neutron stars that emit electromagnetic 
radiation in a beam which, if pointing towards Earth, is 
observed as a "lighthouse" of great regularity. We neglect 
the contribution of recycled pulsars [5]. 

Magnetars [TUHT3"] are neutron stars with a magnetic 
field a few orders of magnitude stronger than usual pul- 
sars. That magnetic field may support large ellipticities 
[14] leading to an enhanced production of gravitational 
radiation. We obtain an upper limit for the background 
produced by the magnetars. 

Gravitars [T5J [TS] are hypothetical neutron stars that 
have a magnetic field weaker than usual pulsars, and lose 
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rotational energy primarily via gravitational radiation. 
There may exist a population of gravitars that cannot yet 
be detected because they emit very little or no electro- 
magnetic radiation. A simulation performed in [17] shows 
that the conditions for neutron stars to be gravitars de- 
scribed in [15] are possible. In this paper we investigate 
the detection prospects for the background produced by 
such a population. The ensemble of gravitars provides an 
upper limit for the background of rotating neutron stars. 

This work is a follow-on study to [TB] , where the back- 
ground produced by binary systems is studied (including 
binaries formed by white dwarfs, neutron stars, and black 
holes). With both papers, two of the most promising 
sources of contemporary background are covered. 

Other potential sources of contemporary background, 
not discussed in this paper or in |18| . are newborn neu- 
tron stars undergoing r-mode instabilities [19), [20] , com- 
pact objects captured by massive black holes [3T], inspi- 
ralling black hole binaries with intermediate or extreme 
mass-ratio |22j . supernovae |23j . and population II and 
III stars [241 ES]. 

Besides the contemporary background, there may exist 
a primordial one [SJ 1261 127) , arising from processes in the 
early history of the universe. 

We calculate what part of the total background of ro- 
tating neutron stars is unresolvable (commonly named 
confusion noise or stochastic background). The signals 
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composing this part cannot be distinguished from each 
other or subtracted from the data of a gravitational wave 
detector (we do not study the problem of the subtrac- 
tion of resolvable signals, treated, for example, in [2"B1 - 
[3U]). The resolvability of the background is quantified 
by the overlap function, Af{f, A/, z), introduced in |18j . 
This function gives the expected number of signals, with 
redshifts smaller than z, that are observed within a fre- 
quency bin [/, / + A/], where A/ is the frequency res- 
olution allowed by the detector and the data analysis 
method. When a frequency bin is constantly occupied by 
one or more overlapping signals, i.e., Af(f, A/, oo) > 1, 
these signals cannot be disentangled, and form an unre- 
solvable background. 

The spectral gravitational wave density parameter, or, 
simply, spectral function, is often used to quantify 

the background [T] . It gives the average energy density of 
gravitational radiation (per logarithmic frequency inter- 
val) divided by the critical density. The generalized spec- 
tral function [TB], A/,A/"o), has the same meaning as 
J7(/), but it quantifies only the part of the background 
with more than Ao overlapping signals per frequency bin. 
The total background and the unresolvable one are cal- 
culated by taking Af — and Af = 1, respectively. In 
this paper, the spectral function accounts for the time 
evolution of the systems, that is not assumed to be short 
compared to cosmic time scales. 

Previous work has studied the gravitational wave back- 
ground from pulsars [3T] and magnetars [52"H51] . These 
articles assume that all neutron stars are formed with 
the same initial spin frequency. We show that the results 
change dramatically if the initial spin frequency follows 
a probability distribution. In particular, for one of the 
distributions considered [35], the detection of the back- 
ground by present and planned detectors is rather unre- 
alistic. 

For some of the models considered, the detection of 
the background of rotating neutron stars could be pos- 
sible by cross-correlating two interferometers of the Ein- 
stein Gravitational Wave Telescope (ET), assuming two 
of the proposed configurations (called ETB and ETD) 
[55] , Furthermore, this background is unresolvable. The 
current generation of present ground-based detectors |37l - 

, and the advanced version of the Laser Interferome- 
ter Gravitational Wave Observatory (aLIGO) [40], are 
not sensitive enough to detect this background. For fu- 
ture space missions like the Big Bang Observer (BBO) 
[2"5] and the Decihertz Interferometer Gravitational Wave 
Observatory (DECIGO) [41], the detection is rather un- 
likely. 

The outline of the paper is as follows: 

In Section [II] the notation and nomenclature of the 
paper is explained, and the quantification of the grav- 
itational wave background, its resolvability and de- 
tectability are briefly reviewed. A general formula for 
0(/, A/, Afo) is derived for a population of systems that 
emit at different times and locations, without assuming 
that the evolution of each system is short compared to 



cosmological time scales. We also give a formula for 
A/o(/, A/, z) which is more general than the one pre- 
sented in the previous work |18j . 

In Section |III| we expand upon the expressions of 
fi(/, A/,A/o) and Ao(/, A/, z), to account for the evo- 
lution of the population. We obtain formulas that de- 
pend on the energy and frequency evolution of a sys- 
tem, the initial frequency distribution and the formation 
rate of the ensemble, and certain cosmological parame- 
ters. Then, assuming that all systems start emitting at 
the same frequency and evolve in short time scales, we 
obtain the formula of the spectral function that is com- 
monly used in the literature. 

In Section IIVI we describe the models assumed for a 
neutron star and its population. 

Section |V| contains the main results of the paper. We 
present a robust upper limit for the background of rotat- 
ing neutron stars, the gravitar limit. We then obtain the 
background produced by gravitars, and study the like- 
lihood of planned detectors to observe it and to place 
limits on the abundance of gravitars. The most realistic 
expectation of the background of rotating neutron stars is 
calculated, using a magnetic field and an ellipticity distri- 
bution from the literature. An upper limit on the back- 
ground produced by magnetars is obtained. We study 
the detection prospects of ETB, ETD, BBO and DE- 
CIGO, for different assumptions on the initial frequency, 
magnetic field, and ellipticity of neutron stars. 



In Section VI we compare our results with others from 
the literature. We also comment on the insensitivity of 
the spectral function on the choice of star formation rate. 

The main results and conclusions are put together in 
Section VII First, in Section VII A| the technical achieve- 
ments regarding the calculation of £l(f, Af, Ao) are sum- 
marized. Then, in |VIIB| we compress all results and 
predictions regarding the detection of the background of 
rotating neutron stars. A non-specialized reader inter- 
ested only in the main conclusions should read the latter 
section. 

In Appendices [X] and [B] we give simple approximate 
formulas for the spectral function of the background of 
rotating neutron stars, and also for the one of binary sys- 
tems. Finally, in Appendix [C] we point out a feature in 
the gravitar limit that is analogous to Blandford's argu- 
ment ITol. 



II. CHARACTERIZATION OF THE 
BACKGROUND: AN OVERVIEW 

We follow the notation and terminology explained in 
Section II of [18] . The index "e" (for emitted) is used 
for frequencies and energies of the gravitational waves, 
as well as intervals of time, measured close to the sys- 
tem (for example, a single rotating neutron star) at the 
time of emission of the radiation. Observed frequencies, 
energies and intervals of time (measured here and now) 
have no index. Emitted quantities f e , E e , and Ai e (and 
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infinitesimal emitted intervals df ei dE e , and dt e ) are af- 
fected by the expansion of the universe. They are related 
to the observed quantities /, E, and At (df, dE, and dt), 
by 

/=[l+z]- 1 /e, df=[l + Z ]- 1 df e , (1) 



and 



E=[l + z]~ 1 E e , dE = [1 + zY x dE e , (2) 



At = [1 + z]At e , dt = [1 + z]dt e , (3) 



where z is the cosmological redshift. Any given function 
x that depends on f e can be written in terms of observed 
frequencies. The notation x| . means that the function 
x(f e ) must be written in terms of observed frequencies, 
i.e. x\ f = x(f[l + z]). 

For convenience, a lookback time interval is sometimes 
used, and denoted by an index L. The relation between a 
lookback time interval At L and an ordinary lookforward 
time interval At is At L — — At. 



A. Quantification of the background 

The gravitational wave background is usually charac- 
terized by the spectral energy density parameter [I] (or, 
simply, spectral function), 



«(/) 



pin(f) £ln(/) 



C 2 p c 



(4) 



where c is the speed of light. The present critical density 
of the universe is 



Pc 



8nG ' 



(5) 



where G is the gravitational constant, and Hq is the 
present Hubble expansion rate, of 74.2 kms^ 1 Mpc^ 1 
[4"2l |4"3"] . The function £i n (/) is defined in such a way 
that £i n (/)c?ln/ is the energy per unit volume of gravi- 
tational waves between In / and In / + d In /. Thus, fl(f) 
is related to the total density of gravitational radiation 
in the universe, that is 



Pgw 



p ia (J)dbxf = p c fi(/)dln/. (6) 



other hand, one loses information when using the spectral 
function for a resolvable background. 

We now derive for an ensemble of many sources, 
emitting at different times and locations, that can expe- 
rience a time evolution. The radiation we observe today 
has been produced by many individual systems in the 
past. The energy emitted by one system during an in- 
finitesimal interval of time is 



dE, 



dE e 
dt e 



-dtp 



dE e dt e 
dt e dti 



dti 



dE e 
dtp 



-dti 



(7) 



Two waves that reach us now and were emitted at differ- 
ent lookback times t\ and t\ + dt% , have different red- 
shifts z and z + dz. Lookback time intervals can thus be 
written as redshift intervals, 



dti 



dti 



-dz . 



(8) 



The number of systems, per unit comoving volume, con- 
tributing to the background with observed frequencies 
between In / and In / + d In / is 



dn 



dn 
dh^f 



din/. 



(9) 



The present energy density of gravitational waves, per 
unit logarithmic frequency interval, produced by the col- 
lection of all systems is 



em(/) = 



*° dE dn 
~dt d\nf~ 



dt = 



dE dn dt L 
dt d In / dz 



dz, (10) 



where to is the current age of the universe. Using Equa- 
tions Q, (|]), and @, we can write 



[1 + z 



dtp 



dn 



: d In f e 



dt^ 
- dz 



dz . 



(11) 



Replacing (11) in Q , we finally reach the formula for the 
spectral function of the total background, 



«(/) 



dE e 
dtp 



dn 



,d\nf e 



■ dz 



dz. (12) 



The functions dE e /dt e and dn/dhif e are obtained in 
Sections |IV C| and |IIIB| respectively, for the ensemble 
of rotating neutron stars. The function dtp /dz depends 
on the choice of the cosmological model; we assume a 
Lambda-Cold Dark Matter universe, so 



Here, f2(/) is the spectral function of all sources of grav- 
itational radiation in our past light cone. For simplicity, 
we use the same symbol to characterize the background 
produced only by the systems we are interested in (ro- 
tating neutron stars). 

The spectral function fully characterizes a Gaussian, 
stationary, isotropic and unpolarized background pp. As 
claimed in [TS], the spectral function is also the right 
tool to characterize an unresolvable background. On the 



where 



dti: = 



[l + z]H £(z) 



dz , 



£(z) = y/n m [l + z} 3 + «A • 



(13) 



(14) 



Here, O m and Oa are the density parameters of matter 
and dark energy, respectively, whose values [44] are as- 
sumed to be f2 m = 0.27 and Ha — 0.73. In Section [II B| 
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we explain how to modify the integration limits of Equa- 
tion (12 1 to measure only the unresolvable part of the 



background. 

In the literature, one usually finds the spectral function 
written as 



i 



dE P 



dln/e 



dt L 

h(z)—^dz 
dz 



(15) 



(see, for example, Equation (35) of [18 or Equation (5) of 
[45]). In Section IIIC we show that (15) can be derived 
from ( |12[ ) if one assumes that systems are short-lived and 
start emitting with the same initial frequency. 



B. Resolvability of the background 

The overlap function, Af(f,Af,z), allows us to define 
and quantify the resolvability of the background. We now 
define the overlap function more generally than in |18j . 



r f f+Af , dn 

JO Jf 1 + Z df e 



, dz 



jdfdz'. (16) 



Here, dV c is the element of comoving volume, given by 



dV r = 4tt 



H S{z') 



dz' 



H S(z) 



dz . 



(17) 



The frequency resolution Af can be chosen as the in- 
verse of the observation time (typically of order one year) . 
The condition of unresolvability is fulfilled from a cer- 
tain redshift z, i.e. Af(f,Af,z) > 1, when each bin is 
always filled by one or more signals. These signals can- 
not be distinguished, because for that we would need to 
improve our frequency resolution; we therefore say that 
they are considered unresolvable |46j . One can invert 
A"(/, Af, z) = Ao with respect to z, obtaining a function 
z = z(f, A/, Ao). Signals with redshifts greater than this 
produce an overlap greater than Ao- Using this redshift 
function as lower limit of the integral in Equation (12), 



one considers only the contribution to the background of 
those signals that produce an overlap greater than Ao. 

Therefore, the spectral function of a background with 
more than Ao signals per frequency bin [/, / + A/] is 



fi(/,A/,Ao) 
1 



, , [1 + 3 

PcC Jz(f,AfJfo) 



_ x dE e 
dt P 



dn 



d In / e 



. dz 



dz. (18) 



For simplicity, we assume that the background is unre- 
solvable if the number of overlapping signals in a fre- 
quency bin is > 1 (other criteria are discussed in Section 
V.D of 18 ). Then, the spectral function of the unresolv- 
able part of the background is given by 



olvabb 



,(/) = «(/, A/, 1). 



(19) 



On the other hand, the spectral function of the resolvable 
part is 

^resolvable \f) — ^total(^) ^unresolvable (,/") ■ (^^) 

Here, the spectral function of the total background 
(which coincides with Equation (12)), is 



n 



total 



(/) = fi(/,A/,0) 



(21) 



where the value of Af becomes irrelevant. 



In Section III C| we prove that Equation ( 16 ) leads to 



the definition of the overlap function given in Equation 
(41) of [18], if one assumes all systems start emitting 
with the same initial frequency and the evolution of each 
system is short compared to cosmic time scales. 

The definition of resolvability given in this section is 
the one that was introduced in [18] and that will be used 
throughout the entire paper. More thorough definitions 
can be a topic for future work, for example taking into 
account the ability of the data analysis method to dis- 
tinguish individual signals from the instrumental noise, 
or to even distinguish unresolvable signals with differ- 
ent amplitudes or directions of arrival. The advantage 
of our definition is that the resolvability becomes an in- 
trinsic property of the background, i.e., independent of 
the characteristics of the detector (such as its sensitivity) 
and of the data analysis method. The only parameter re- 
lated to the observation that affects the resolvability is 
the frequency resolution. However, the observation time 
T bs of any realistic experiment is of the order of a year 
or a few years; the best frequency resolution achievable 
(calculated as Af = T~q), can thus be considered equal 
for all possible detectors. 



C. Detectability of the background 

In practice, the instrumental noise of a detector can- 
not be modeled with perfect accuracy; if an unresolvable 
background is present in the data of a detector, it is there- 
fore indistinguishable from instrumental noise (unless one 
can construct a null stream, which is a very advantageous 
feature of ET 07] )• The usual technique to cope with this 
issue is the cross-correlation of the data of two detectors 
(see, for example, Section 7.8.3 of 05]). 

If a background (characterized by a spectral function 
f2(/)) is present in the data of two interferometers, one 
can cross-correlate the outputs of both, that span an in- 
terval of time T bs- Doing this, the resulting signal-to- 
noise ratio (Equation (7.241) of 05]), or SNR, is given 
by 



SNR=^f 

47T 2 



2T Q 



bs 



df 



-,1/2 



(22) 



Here, S n ,i(/) and S ny 2(f) are the noise spectral densities 
of the detectors 1 and 2, respectively, and T(f) is the non- 
normalized overlap reduction function (Equation (7.226) 
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of gS]), defined by 



r(/) = ^ /„ d0 L sm{9)ddcos c 



JO 

2tt 



/ # £ F?(9,cf>,i>)F$(9,<l),TP). (23) 



In this definition, Aa; = — xi, where is the position 
of the detector d, and u{9, <p) is a unit vector pointing to 
the direction defined by the angles 9 and </>. The function 
F%(9,cf>,ip) (for the detector d =1 or 2, and for the po- 
larization p — + or x ) is the antenna pattern function, 
evaluated at the direction (9, 0), for a wave with a polar- 
ization angle ip. The antenna pattern functions can be 
found in Section II.B of [47] for ET, and in Section 2.1 of 
[4Ti] for aLIGO. Notice that Equation (22 1 is equivalent 
to Equation (3.75) of PQ; however, the normalized over- 
lap reduction function j(f ) defined in [48] and in [1] are 
only equivalent for the case of an L-shaped detector. A 
detailed study on the overlap reduction function can be 
found in 1501. 



Following [48 , the non- normalized overlap reduction 
function can be written as 



r(/) = Fi, 27 (/) . 



(24) 



For two colocated and coaligned detectors, 7(7) = 1 
for all frequencies. For the correlation between two in- 
terferometric V-shaped detectors like ET, one obtains 
•Pi, 2 = 3/10, whereas for L-shaped detectors like aLIGO, 
-Pi, 2 = 2/5. In Section [V] the SNR is calculated for two 
interferometers of ET, and for two aLIGO interferome- 
ters (one at Livingston and one at Hanford), using the 
full overlap reduction function (Equation (23)) and as- 
suming an observation time of one year. For simplicity, 
the SNR for two interferometers of BBO or DECIGO will 
be calculated by using Fi,2 = 3/10 and 7(/) = 1, and an 
observation time of one year. 

Figure hi shows the spectral strain sensitivity \J S n (f) 
of aLIG(J [5T], two possible configurations of ET [52"] . 
DECIGO [53 and BBO [S3]. We consider two pairs of 
detectors: the two aLIGO detectors, at Hanford and Liv- 
ingston, and two V-shaped ET detectors sharing one arm 
of the triangle. 

A background is said to be detectable if it produces 
SNR larger than a certain threshold value. Be aware 
that a background can be resolvable without being de- 
tectable; it would consist of signals that are separated in 
a frequency-time plot but would be buried in noise (for 
example, instrumental noise, or confusion noise produced 
by another background). 



-22 



-23 



-24 



CD 
O 



\\ DECIGO 



aLIGO 



BBO 



ETTV 



ET-B 



-25 

-10 12 3 

log 10 (//Hz) 

FIG. 1. Spectral strain sensitivity of aLIGO, two possible 
configurations of ET (named ETB and ETD), DECIGO and 
BBO. The sources of the curves are given in the text. 



III. CHARACTERIZATION OF THE 
BACKGROUND: A DETAILED DERIVATION OF 

n(f,Af,Af ) FOR AN EVOLVING POPULATION 
OF SYSTEMS 

A. Formation rate of systems 

The comoving density rate of systems formed h(z) (or, 
simply, rate), is defined such that h(z)dz is the number 
of systems formed per unit emitted interval of time, dt e , 
per unit comoving volume, dV c , between redshifts z and 
z + dz. 

Sometimes it is convenient to write the rate as a func- 
tion of time, instead of redshift. We define a function, 
T(z), that gives the interval of time elapsed between the 
formation of the first systems (at redshift z max ), and the 
formation of the systems at redshift z. This function can 
be derived, for our cosmological model, using the formu- 
las given in Section 11.13 of [54] , 



T(z) 



asinh 




Here, we have introduced the constant 



£ = asinh 




[1 ~~~ ^max] 



-3/2 



(25) 



(26) 



that imposes a time offset between the Big Bang and the 
formation of the first systems [SS] ■ One can invert T{z) = 
At with respect to the redshift and obtain another useful 
formula, 



Z(At) = 



Sinn 

n A I 




-1 -2/3 



- 1 , 



(27) 

This gives the redshift observed in a signal that was emit- 
ted an interval of time At after the formation of the first 
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systems. Using Equation (|27j), one can write the rate as 
a function of time, n(Z(t)). 



at its formation is given by f- m i = f- m i(f e ,t e ). Using this 
function we can rewrite Equation ( 29 1 , 



Time evolution of the ensemble 



9/ini 

Pe{fe,t e ) = Pini (/ini (fe, *e)) -jjj- (fe, *e) . 



(30) 



We now 
[dn/dbif, 



explain how to calculate the term Equations ( 28 ) and ( 30 1 can be combined into 



/ 



in Equation ( 12 ), that is the number of sys- 



tems per unit comoving volume per unit logarithmic fre- 
quency interval emitting around an observed frequency 
/• 

Suppose a gravitational wave of redshift z was emitted 
by a system an interval of time T(z) after the formation 
of the first syste ms (re call the definition of the function 
T(z), in Section III A). At the instant of emission, the 



system had already evolved during a certain interval of 
time t e (smaller than T(z)). The system was thus formed 
an interval of time T(z) — t e after the formation of the 
first systems. At the instant of formation, the system was 
emitting waves that have now a redshift Z(T(z) — t e ). 
The formation rate at that instant was n(Z(T(z) — t e )). 
Then, during an infinitesimal interval of time dt e , the 
number of systems formed per unit comoving volume that 
emit waves of redshift z is h(Z(T(z) — t e ))dt e . 

Assume a probability density function p e (f e ,t e ), such 
that p e (fe,te)dfe is the probability of a system to emit 
between f e and f e + df e after a time evolution t e . Then, 
the number of systems formed during dt e per unit comov- 
ing volume that emit waves of redshift z in the frequency 
interval [f e ,f e + df e ] is p e (f e ,t e )h(Z(T(z) - t e ))dt e df e . 
The total number of systems per unit comoving volume 
per unit frequency interval emitting with redshift z and 
frequency f e is 



dn 

We 



T(z) 



Pe(fe,t e )n(Z(T(z) -t e ))dt e 



(28) 



We now show how to calculate the probability density 
function p e (f e ,t e ). For this purpose, we follow a similar 
approach as [161 156] , although we do not use any distri- 
bution of systems in the galaxy, but rather assume that 
systems are homogeneously distributed in the universe. 

Additionally, suppose we know the initial probability 
density function p in i(/ini), such that Pini(/mi)d/ini is the 
probability of a system to emit between /; n j and fi n i+df inl 
at the instant of formation. Now we make the following 
assumption: a formed system never stops emitting grav- 
itational waves. This means that all systems that were 
initially emitting in the frequency range [/;„;, /i n ; + df m i] 
are now emitting in [f e , f e + df e ]. Then, in order to con- 
serve the number of systems, 



Pe(fe,t e )dfe = Pini (/ini) df u 



(29) 



must be fulfilled. 

The radiation we observe now from a system, at fre- 
quency /, was emitted in the past at frequency / e ; that 
system was formed an interval of time t e before emitting 
at f e . The frequency at which the system was emitting 



dn 

We 



T{z) 



dt e pini (/ini {fe > ^e)) 



x d -^{f e ,t e )n{Z(T{z)-t e )). 



(31) 



Finally, we can rewrite Equation (31) in terms of ob- 
served frequencies per logarithmic frequency interval to 
obtain 



dn 





. dn 


/ 


dfe_ 



f 

--f[l + z] / dt eP - ml (f ini (f[l + z],t e )) 
Jo 



dfe 



'(f[l + z],t e )h(Z(T(z)-t e )). (32) 



In Section |III C[ we derive general formulas for the spec- 
tral function and for the overlap function, that depend on 
the shape of the function Pini(/ini); then, we highlight a 
special case where all systems emit at the same frequency 
at the instant of formation. 



C. Spectral function and overlap function 

1. General case 

The spectral function of a background, with more than 
Ao signals per frequency bin A/, produced by an ensem- 
ble that follows an initial frequency distribution J3i n i(/ini)j 
is given by 



n(/,A/,M)) 



dz- 



dEe 

dtp 



d^ 
dz 



PcC* Jz{fAfMo) uu e f 
fT(z) 

x / dt ePini (f ini (f[l + z],t e )) 
Jo 

x^(f[l + z],t e )h(Z(T(z)-t e )). 



(33) 



This is obtained by replacing Equation (32 1 in (18); the 
upper limit of the integral has been replaced by z m ax, 
since we assume that no systems were formed at larger 
rcdshifts. 

The overlap function of a background produced by an 
ensemble that follows an initial frequency distribution 
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PiniCfini) is given by 



rz rf+Af ,y 

Af(f,Af,z)=J dz'J d f'[l + z'}^ 

rT(z') 

x / dt e p- mi {f ini {f[l + z%t e )) 
Jo 

x ^(f[l + z>],t e )h(Z(T(z')-t e )). 

(34) 

This is obtained by replacing [dn/df e ]\f (from Equation 



(32 1) in (16 1. The frequency resolution is typically much 



smaller than the range of frequencies of interest for rotat- 
ing neutron stars. Equation (34) can thus be simplified 
with the approximation 



JV(/,A/,z)« f dz'Af[l 
Jo 



j,dK 
dz' 



z 



/■TO') 

x / dt ePini (f ini (f[l + z'],t e )) 
Jo 

x d ^(f[l + z') ) t e )n(Z(T(z')-t e )), 

(35) 



which is accurate as long as Af <C /. 



2. For a fixed initial frequency 

From all possible initial frequency distributions 
Pini(/ini), we now study a particular case with the form 



Pini(/ini) = <K/im - /fix) 



(36) 



Using this distribution, one assumes that all systems 
start emitting at a fixed initial frequency fa x . Suppose 
we have a function /i n i(/e,0 (that was introduced in 
Section IIIIBl and will be derived in Section IIVDI for the 



case of a rotating neutron star) . Using one of the proper- 
ties of the Dirac delta function, one can rewrite Equation 



(36) as 



Pini (/ini (/e > ^e)) 



dfir 



(/e,T e (/ nx ,/ e ))<5(i e - r e (/ fix ,/ e )). 



Here, T e (/ nx ,/ e ) is the interval of time that a system 
spends emitting between / nx and f e such that 



/ini {fe ) 7~e (/fix 5 fe) ) /fix 



(38) 



is fulfilled. The transformation of the Dirac delta func- 
tion performed is valid as long as df ln i/dt e is non-zero for 
all values of f e and t e . This condition holds for rotating 
neutron stars [57|. By replacing Equation (37) in (31), 
one obtains 



Qt e d /i n i 

" "(/e, Te(/fix,/e))^7^(/e, Te(/fix, fe)) 



dn 

dfe <?/ini Of e 

x n(Z(T(z) - r e (/ fix , f e )))B(z, f e ) , (39) 



where we have introduced the function 

Q(z, f e ) = 8(T(z) - r e (/ fix , / e ))0(r e (/ fix , / e ) - 0). (40) 

Here, 9(x — y) is the Heaviside step function; it is equal to 
one for x > y, and equal to zero for x < y. For rotating 
neutron stars, dt e /df e is not a function of t e . Therefore, 
without loss of generality, we write dt e /df e (f e ). Then, 



dn dt e 
dje = dfc 



(f e )h(Z(T(z) - T e (/fix, fe))Mz, fe) . (41) 



We now substitute this result in the formulas of the over- 
lap function and the spectral function. 



Inserting Equation (41) in (16), we obtain 



r f+Af dt 

dz' / d /'[l + z']^(/'[l + z']) 

Jf dfe 



dVr 



x -^n(Z(T(z') - r e (/ fix , f[l + z'])))Q(z>, f[l + z'}) , 

(42) 



which is the overlap function of a population of systems 
that start emitting at the same initial frequency / nx . 
Similarly, inserting Equation (41) in (18), we obtain the 
spectral function 



fi(/,A/,M)) 



1 



PcC Jz(f,Af,M,) 



dE P 



dlnf e 
dtl 



xh(Z(T(z) - r e (/ fix , f[l + z])))Q(z, f[l + z\)^dz . 



(43) 



Performing the same approximation as in Equation ( 35 1 , 
the overlap function can be simplified as 

Af(f, A/, z) « Af£ dz'[l + z']^^(f[l + A) 



x h(Z(T(z') - r e (/ nx , f[l + z'])))6(z', /[l + z']) 



(44) 



(37) which is accurate for A/ <C /. 



3. For short-lived systems with fixed initial frequency 

We now prove that if systems are assumed to evolve 
rapidly compared to cosmic time scales, then one obtains 
the definition of J\f(f,Af,z) and f2(/, A/, Af ) given in 
Equations (41) and (44) of [13], respectively. Under this 
assumption, T e (/ nx ,/ e ), which is the time a system has 
evolved since formation, is much smaller than T(z), and 
therefore, 

h(Z(T(z) - r e (/ fix , /„))) « n(Z("7»)) = n(z) . (45) 



Using this, we can rewrite Equation (42) to obtain 

,dV c 



A^(/,A/,z) 



dz'n(z') 



dz' 



»[/+A/][l+ Z '] 

< / ^ {fe) e(z'j e )df e 

>f[l+z'] a Je 



*low(/) 



r e (/,A/,z>(z')^^'. (46) 



Here, the function T e (/, A/, z) gives the interval of time 
that a system, whose radiation is now observed with red- 
shift z, spends emitting between observed frequencies / 
and / + A/. Alternatively, inserting Equation (45) in 
( 43 1 , the spectral function becomes 



1 



PcC Jz(f,AfJT ) 



dE 

[1 + z]- 1 dhe 



dkif e 



dt L 
h(z)—^dz . 
dz 



(47) 



In Equations (46) and (47), the functions z\ ovi (f) and 



■*upp 



(/) ensure that the integration is performed only 
where 6(2, f[l + z]) is non-zero. 

Let us examine the spectral function of the total back- 
ground, i.e., the one obtained by imposing no restriction 
(Mo = 0) in the number of signals per frequency bin. By 
inverting Af(f,Af,z) — (Equation (46)) with respect 
to the redshift, one obtains z(f, A/, ffj = z\ ow (f). Re- 
placing this in Equation ( |47[ ), the canonical formula for 
the spectral function (used for example in [T51 HHl EB] ) is 
recovered, 



1 



,(/) 



dE, 



d\nf e 



dt L 
h(z)—^dz . 
dz 



(48) 



By carefully studying the limits of this integral, one re- 
alizes that the redshift functions z\ ow (f) and z upp (f) are 
not exactly the same as the ones defined in Equations 
(37) and (38) of [18], or in Equations (10) and (9) of [58]. 
The difference, however, occurs only at the low-frequency 
part of the spectrum, at which the time scales needed for 
the systems to evolve are comparable to cosmic ones [59] . 



IV. MODEL FOR THE ENSEMBLE OF 
ROTATING NEUTRON STARS 

A. Neutron star model 

A neutron star is modeled as a rigid rotating ellipsoid 
of mass to. Its semiaxes with respect to the coordinate 
axes x, y, and z have lengths a, 6, and c, respectively. The 
rotation occurs around the z-axis at an angular velocity 
to, which slowly decreases in time. Assuming a uniform 
density, the moment of inertia about the z-axis is 



(49) 



The ellipticity is defined by 



e=[h-h]/h 



(50) 



where I\ and I2 are the moments of inertia about the x 
and y axes, respectively. The average ellipticity of the 
ensemble of neutron stars is very uncertain; reasonable 
values for e can range from 10 -8 to 10~ 4 . For such val- 
ues, e ss [b — a]/a, and, replacing it in (49), we can very 
accurately approximate 



2m a 2 



(51) 



We assume m = 1.4 M Q (where M is the solar mass) 
and a = 12 km [601 161 j . obtaining a moment of inertia of 
/ w 1.6 x 10 38 kgm 2 (see Sections 3.1.1 and 3.1.2 of [8 
and references therein for a discussion about these val- 
ues). Other mechanisms of gravitational wave emission, 
like dynamical bar-mode [62] or r-mode [19] instabilities, 
are not taken into account in this work. 

A neutron star behaves like a rotating dipole magnet. 
The value of the magnetic field at the magnetic pole is 
B, where it forms an angle a (assumed, for simplicity, 
of a = 7r/2) with the direction of the rotation. We do 
not consider any magnetic field decay [53] ; B is the same 
during the entire life of the star. 

With this toy model, one can infer the average value of 
B by looking at the rotating period and its time deriva- 
tive (the so-called P — P diagram) of a pulsar catalogue 
[6"1] . Without taking into account recycled pulsars, a 
reasonable average value for pulsars is B = 10 8 T. For 
magnetars, larger values (of B pa 10 10 T) can be reached. 

The maximum gravitational wave frequency at which 
a rotating neutron star can emit is estimated by 



f 

J ] 



escape 

max 



1 Gm 



3.3 kHz. 



(52) 



Above this frequency, the material at the equator would 
have enough velocity to escape the gravitational poten- 
tial, since the latter becomes lower than the centrifu- 
gal potential [60]. A gravitational wave frequency of 
3.3 kHz corresponds to a rotation period of 0.6 ms, which 
is roughly the half of the fastest rotation period known 
in a pulsar [65]. A more realistic estimate [66] of the 
maximum frequency is 



pRochc 
J max 



-1 3/2 



f 

■I 1 



escape 
max 



1.8 kHz. 



(53) 



This frequency takes into account the deformation of the 
equatorial radius because of the rotation. From now on, 
wc make the choice 



fn 



escape 
max 



(54) 



which leads to the most optimistic results, regarding the 
detection of the background. The conclusions of this pa- 
per would be unaffected, however, by choosing the alter- 
native maximum frequency in Equation ( 53 1 . 
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B. Formation rate of neutron stars 

The amount of mass converted into stars per unit emit- 
ted interval of time per unit comoving volume between 
redshifts z and z + dz is given by p(z)dz, where p(z) is 
the star formation rate. The models for the star forma- 
tion rate [67-74 usually present a similar shape: p(z) 
increases from its local value (at z = 0) until z « 1 or 2, 
and then decays, reaching negligible values for redshifts 
larger than 5 or 6. For this reason, the range of redshifts 
considered in the calculations is [0, 2 max ], with 

z max = 5 . (55) 

All calculations shown in the plots of Section |V| are ob- 
tained by assuming the star formation rate given in Sec- 
tion 5.4 of [Ml, 



bz 



l+[f] 



1 Mgyr^Mpc" 



with the parameters found in Table 
namely (a,b,c,d) = (0.0170,0.13,3.3,5.31 
iVflOOkms^Mpc- 1 ] = 0.742. In Section[VIB 



I of 

and 



(56) 



ment on the fact that the specific choice of star formation 
rate does not affect the spectral function significantly. 
Furthermore, the results would not be affected by the 
use of a constant rate. 

Only a fraction A of all stars formed become neutron 
stars, so the rate is 



h(z) = Xp(z) . 



(57) 



The fraction of stars formed with masses between m and 
m+dm is $(m)<im, where <&(to) is the initial mass func- 
tion. We assume a Salpeter initial mass function |75) . 

Am- 2 35 , 



$(to) 



(58) 



where the value of the normalization constant A turns 
out to be unnecessary, as we now see. The fraction of 
stellar mass converted into neutron stars is 



A = 



r20 Mq 

J8M 



$(m)dm 



Jo.iAfc "»#(m)dm 



= 5 x 10" 3 M, 







(59) 



The denominator is the average mass of a star (consider- 
ing all stars with masses between 0.1 M and 100 Mq), 
and the numerator is the fraction of stars that can be 
progenitors of neutron stars (namely stars with masses 
between 8 Mq and 20 Mq). The value of A tells us that, 
for each 10 3 M of gas converted into stellar mass, 5 neu- 
tron stars are produced. We assume that this fraction is 
the same at all redshifts. 



C. Energy evolution 

The rotational energy of a system is given by 

1 



E,, 



(60) 



For convenience, we write the angular velocity ui e in 
terms of the frequency of the emitted gravitational waves, 
/ e , which fulfills 



7I"/e 



(61) 



(61 1, one obtains 



Differentiating (60) with respect to the time, and using 

dfe 



dt e 



dt f 



(62) 



In the following we show that df e /dt e is negative, thus, 
dE erot /dt e is also negative; the system loses rotational 
energy with the time. This energy loss is due to the emis- 
sion of both electromagnetic and gravitational radiation. 

Suppose that the system lost energy only via the mag- 
netic dipole emission. Such a system is studied in 0]. 
Rewriting Equation (15) of that paper with our notation, 



dE e , Bld _ 2n 5 a e B 2 sm 2 {a) 
dt e 3c 3 /io 



ft = -^ft, (63) 



where the index 'md' stands for magnetic dipole and 

3^ c 3 / 



3md 



4n 3 B 2 a 6 sin 2 (a) 



(64) 



Here, po is the magnetic permeability of the vacuum. 
We have used that B = poH, where H is called R\ in 
[4]. The magnetic field is often expressed in Gaussian 

units [33l El] . For clarity, recall that -Bcauss = 



B. 



Equation (63) gives the rate at which the system loses 



energy by emitting electromagnetic waves. 

Suppose now that the system emitted only gravita- 
tional radiation. This system is studied, for example, in 
Section 4.2.1 of @8]. Rewriting Equation (4.227) of [4"8] 
with our notation, 



dE e 



dt f 



32tt 6 G€ 2 I 2 
5c 5 



ft 



Try 

"4A 



-ft 



(65) 



where the index 'gr' stands for gravitational radiation, 
and 



5c 5 



1287r 4 Ge 2 r 



(66) 



Equation ( 65 1 gives the rate at which the system loses 



energy by emitting gravitational waves. 

The system we study loses energy via both magnetic 
dipole emission and gravitational radiation. Thus, the 
total loss of energy (which is a loss in rotational energy) 
fulfills 



dE P 



dE e 



dE P 



dt P 



dtf 



dtf,. 



(67) 



Rewriting Equation (67) in terms of Equations ( 62 ) , ( 63 1 , 
and (65), we obtain 



dfe 

dt r< 



1 



25, 



■I e 



md 



i 

4<L 



,5 

J e ' 



(68) 
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where we see that the frequency of the rotation (written 
in terms of gravitational wave frequencies) decreases with 
time. The energy lost by the rotating star is the energy 
gained by the gravitational waves, so 



dEe 

dt„ 



dE„ 



dt P 



HI 1 

A5, 



-ft 



(69) 



is the amount of gravitational wave energy produced by 
a system per unit time. The amount of gravitational 
wave energy produced by one system contained in an 
infinitesimal logarithmic interval of frequency d In f e is 



dE e dEe 
dlnfe ' e df e 



= /e 



dE e 


\dfe~ 


-l 


dt e 


dt e _ 





(70) 



where we have used the chain rule. Using Equations (68) 
and ((69b, 



dEp 



dlnfe Pc + fV 
where we have introduced the cut frequency. 



fc 



26^ 

Smd 



(71) 



(72) 



At this frequency, both terms on the right side of Equa- 
tion (|68|) become equal; this is the frequency at which 



both mechanisms of energy loss "cut" each other. The 



absolute value in Equation ( 70 ) is used because dE e must 



be a positive quantity; it represents the amount of grav- 
itational wave energy within a logarithmic frequency in- 
terval, regardless of whether the energy of the system 
increases or decreases with the frequency. 

The three main expressions of this section are in Equa- 
tions ( 68 ) , ( 69 1 , and ( 71 1 . They can be rewritten, using 
0, in terms of observed frequencies. The resulting for- 
mulas are: 



dfe 

dt R 



—^/ 3 [l + *] 3 --^/ 6 [l + *]' 
2d m d 4dgr 



and 



dEe 

dt e 
dEe 



TT 2 I 



f[± + z] 



dlnfe 



f 



tt 2 I / 4 [1 + z} 4 
f 2 [l + z] 



fl 



(73) 



(74) 



(75) 



respectively. 

We can distinguish two frequency intervals: one where 
the magnetic dipole emission dominates (let us call it 
the md-range) and one where the gravitational radiation 
dominates (the gr-range). The frequency at which both 
mechanisms are equally dominant is fc- For simplicity, 
some of the next calculations are performed in the two 
frequency intervals separately. The energy spectrum can 
be approximated by 



dEe 



dlnfe 



Jn 



(76) 



in the md-range, and by 



dEe 



dlnfe 



ir 2 If 2 [l + zY 



(77) 



in the gr-rangc. 

The braking index rib is defined by [5] 



dfe 
dt e 



-K{feT b , 



(78) 



where if is a constant. Equation (68) shows that the 
braking index is equal to 3 in the md-range and equal 
to 5 in the gr-range for all neutron stars. Observational 
measurements of the braking index, however, obtain very 
different values. For example, in Table 4 of [35] rib is 
smaller than 3 for some known pulsars. Alternatively, 
the braking index measured in other pulsars can be orders 
of magnitude larger than 3, or even negative [76] . The 
results of this paper would differ considerably if one used 
models with different braking indices. The consideration 
of such other models is out of the purposes of this work. 



D. Frequency evolution 

We now calculate the lapse of time r e (/ e ,i, fe,2) that a 
system spends emitting within a certain frequency inter- 
val [/e,ij/e,2]- This is achieved by integrating Equation 
( 68 1 , which leads to the analytical formula 



T~e(fe,l, fe,2) 



9md 



fe.2 fc 



.i+fc 2 ^ 



fe,%Afc + fe,l] 

flAPc + fh] 



(79) 



A system that starts emitting at an initial frequency 
/i„i, needs an interval of time r e (/i n j,/ e ) to reach the 
frequency f e . 

It is useful to obtain a function f lnl = f- m i(f e , At e ), 
that gives the frequency at which a system, that now 
emits at f e , was emitting an interval of time At e be- 
fore. One cannot invert Equation (79) with respect to 



/i n i analytically. For this reason, it is more convenient to 
approximate r e (/ ini ,/ e ) by 



7e(/ini, fe) 

6, 



f-4, _ f-4i 
~ 4 m \ 
^gr[/c _ /ini J 
^md[/ e — /; n i 



+ <W/ e 



fc < fe < fix 
fc 2 } fe< fc < fix 
fe < /ini < /( 



C 



(80) 



With this approximation, one can analytically invert 
T(/ini,/ e ) = At e with respect to / ini , obtaining 
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/ ini (/ e ,Ai e ) 



fc*W. 

fe 2 



-1/4 



A/,. 



At, 



-1/2 



In this equation, the conditions / e < 5/(7 2 + y- 
-1/4 



and f e < 
values for f\ 



At,, 



-1/4 



E. Initial frequency distribution 



/c</e 
1/2 



A/, 



-1/4 



-1/2 



are introduced to avoid unphysical 



fe< 



fc 2 + 



<fe<fc & fe< 
1/2 



2/c 



Ate 



-1/2 



(81) 



I 

[751 [SO]. Another possible distribution, used in [5TL could 
be a Gaussian distribution like that of Equation ( 85 1 with 
fi =50 ms and a =50 ms. This distribution leads to inter- 
mediate results between those of Distributions 1 and 2. 
Regarding gravitars, our current knowledge about their 
population statistics is so poor that any of the previous 
distributions is equally plausible. 



We consider three simple initial frequency distributions 
Pini(/ini) in the calculations. 



The first one was already introduced in Equation ( 36 1 . 
Let us call it Distribution 0. 

Distribution 1 is obtained from the log-normal initial 
period distribution given in |77j . 



Pp(Pini) 



1 



'2iro-P u 



exp 



[ln(iWs) - 
2a 2 



(82) 



where fi = ln(0.005), a = 0.3/log 10 (e) « 0.69. This 
distribution hence assumes that the average initial spin 
period is of 5 ms. The initial period, P ln i, is related to 
the initial spin frequency f- lnl (in terms of gravitational 
wave frequencies) by 



Therefore, 



Pi, 



Pini ( fixxi) 



W pp V f- 

Jin \JB 



(83) 



(84) 



is the corresponding probability density function of the 
initial frequency. 

Lastly, Distribution 2 is obtained from the normal ini- 
tial period distribution given in [35], 



Pp(Po) 



1 



27TCT 



exp 



[Po - A 

2fJ 2 



(85) 



with [i =300 ms and a =150 ms. Similar distributions 
to this one are used in [78] and in [79]. To obtain the 
corresponding probability density function of the initial 
frequency, one can again use Equation ( 84 1 . 



Initial frequency distributions like Distributions and 
1 are more favorable for the detection of the background 
than Distribution 2. Some studies do predict large ini- 
tial frequencies for the population of magnetars [10] ; on 
the other hand, an initial period of 5 ms (like the aver- 
age of Distribution 1) or shorter may be considered too 
small to properly describe the ensemble of known pulsars 



F. Magnetic field and ellipticity distributions 

Some of the calculations in Section [V] are performed 
using a magnetic field distribution and an ellipticity dis- 
tribution. The formulas for the overlap function (Equa- 
tion (35)) and the spectral function (Equation (33)) can 
be modified to take into account these distributions. The 
overlap function becomes 

N(f,Af,z) = 



dep. (e) / dBp B (B)Af(f, A/, z) . (86) 

iB mia 



Inverting Af(f, A/, z) = Afo with respect to the redshift, 
one obtains a function ?(/, A/, A/o). This function can 
be used as a lower limit of the redshift integral in ( 33 1 , to 



obtain £l'(f, A/, No). The spectral function is obtained 
by solving 

A/,A/- ) 



de Pe (e) / dBp B (B)tf(f,Af,tf ). (87) 



In these two formulas, the magnetic field and ellipticity 
distributions are assumed to be independent. We point 
out that these distributions could in fact be correlated; 



as an example, in Section V C we mention that a high 
magnetic field can increase the ellipticity. We now specify 
the magnetic field and ellipticity distributions used. 

The probability density function of the magnetic field 
is taken from [77] (also used in [i~5l ITS] ). 



Pb(B) = 



V2tto-B 



exp 



Kg/T) - A 
2ct 2 



with n = ln(10 8 35 ) and a = 0.4/log 10 (e) « 0.9. This 
means that the average magnetic field is of 10 8 35 T. The 
distribution is normalized to unity between a minimum 
value of -B m in = 10 7,2 T and a maximum value of B max — 
10 98 T. 
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The probability density function of the ellipticity is 
taken from 15J, 



Pe(e) = A- 



exp(-~, 



r[l-exp(-^-)] 
where t is the solution of 



exp(^) -1 



(89) 



(90) 



The values for e and e max are 10 -7 and 2.5 x 10~ 6 , re- 
spectively. The normalization constant A is obtained by 
imposing 



Pe(e)de = l, (91) 
where the minimum ellipticity is e m j n ss 0. 

G. Minimum and maximum frequencies 



In Section |III C 2[ a formula for the spectral function is 
obtained (in Equation ( 43 ) ) , assuming a fixed initial fre- 
quency; Af,M ) contains the function Q(z, f[l+z]), 
defined in Equation (40), that determines the redshifts 



and observed frequencies of the systems that can con- 
tribute to the background. Introducing ([I]) in (40), one 
gets 



e(zj[i + z}) 

= 0(T(z)-T(f Bx ,f[l 



#M/fo/[i + *-o). 

(92) 



We now study the limits that this function sets on the 
possible observed frequencies and redshifts of the grav- 
itational waves, for the ensemble of rotating neutron 
stars. For that, the fixed initial frequency can be re- 
placed by the maximum frequency (in Equation ( |54[ )), 

i*^* /fix /max- 



The first Heaviside step function in ( 92 ) becomes zero 



for a certain observed frequency / = fiow(z). This func- 
tion gives the minimum observed frequency that a gravi- 
tational wave with redshift z can have. Using Equations 
(p5j) and Q, the condition T{z) = r(/ max , f\ ow (z)[l+z\) 



leads to 



, h{z) fc<fx(z) 
flow{z) = <j h{z) h{z) <fc<f, 

f?,{z) /max < fc 



where 



h{z) 



T(z) 



1-1/4 



+ /, 



Mz) = 



<5md 



-4 
max 



f-4 1 
/ maxJ 



[! + *]" 



(93) 



(94) 



fc 2 



-1/2 



[l + z] 



(95) 



and 



T(z) 



Jmd 



+ /, 



-2 
max 



-1/2 



[1 + zY 



(96) 



One should notice that f\ ov/ (z) is an observed frequency, 
unlike fc and / max , that are emitted frequencies; the e- 
index in the two latter quantities has been omitted to 
ease the notation. 



The second Heaviside step function in Equation (92) 



becomes zero when evaluated at the observed frequency 
/ = fupp(z). This function gives the maximum ob- 
served frequency that a gravitational wave with red- 
shift z can have. Using Equation (80), the condition 
T(/max, f U pp(z)[l + z] ) = leads to 



/upp(^) — /max[l "i" z] 



(97) 



Again, notice that / upp is an observed frequency, whereas 
/max is an emitted frequency. 

With the previous results, we can calculate the max- 
imum and minimum observed frequencies possible. The 



maximum observed frequency, as Equation ( 97 1 clearly 



shows, is achieved at redshift 0, and is precisely / ma x- 
On the other hand, to find the minimum observed fre- 



quency, one has to minimize Equation ( 93 1 . The redshift 
at which f\ ow (z) is minimum is the solution of 



nN <5md 
z ) + 72— 



P 

J 17 

if fc < h(z), of 



dz 



(z) + 2[l + z]- 1 =0, (98) 



T(z) + -^ + 



4<5gr /, 



4 

max 



dT 

dz 



(z) + 2[l + z]- 1 =0, (99) 



if fi(z) < fc < /max, and of 



T{z) + 

J max 



-i -1 



dT 

dz 



(z)+4[l + z]- J =0, (100) 



if / max < fc- As a good approximation, one can assume 
that T(z) > <S md / m ax_and T{z) > S ST f~^. Doing this, 
Equations (98) and (100) depend only on cosmological 



parameters, and their numerical solutions are 



and 



3.39, 



z md » 1-54, 



(101) 



(102) 



respectively. The solution of Equation ( 99 1 will depend 
on the values of the astrophysical parameters (-B, e, et 
cetera), but must lie between z m( j and z gr . As an ex- 
ample, for a rotating neutron star with B = 10 8 T and 
e = 10~ 7 , the cut frequency is fc > /max, so the mini- 
mum observed frequency is given by /3(z m d) ~ 86mHz. 

The upper plot in Figure [2] shows the redshifts and ob- 
served frequencies that the gravitational waves from the 
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FIG. 2. Redshift versus observed frequency of the gravita- 
tional waves produced by the ensemble of gravitars, assuming 
an ellipticity of e = 10~ 7 . The vertical axis gives the red- 
shift of the gravitational waves observed today (redshift at 
present). Those waves were emitted by gravitars that started 
radiating at the redshifts given by the gray scale (redshift of 
formation). The solid lines follow points of equal redshifts of 
formation (corresponding, from bottom to top, to redshifts 1, 
2, 3, 4, and 5). The upper plot accounts for all signals in the 
universe, whereas the lower plot accounts only for unresolv- 
able signals. 
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FIG. 3. Same as Figure |2j but for the ensemble of rotating 
neutron stars with magnetic field B = 10 8 T and ellipticity 
e = 10~ r . 



ensemble of gravitars with ellipticity e = 10 -7 can have. 
The lower plot is the same, but only for the unresolvable 
waves (those that produce an overlap larger than Nq = 1 
in a frequency bin A/ = lyr -1 ). A point in the graph 
tells the present redshift observed in a gravitational wave 
emitted by a gravitar, and the gray scale tells the redshift 
of the waves emitted by that same gravitar at the instant 
of its formation. Figure [3] is analogous to [2j but for the 
ensemble of rotating neutron stars with B = 10 8 T and 

e = lfr 7 . 



V. RESULTS 

A. Gravitars 

Let us first give a definition of gravitar: it is a ro- 
tating neutron star that emits gravitational waves at a 
frequency f e > fc, i.e., the dominating mechanism for 
the loss of rotational energy is the emission of gravita- 
tional waves (see the definition of the cut frequency in 
Section liVCl). 



The gravitar limit is an upper limit on the gravita- 
tional wave background produced by rotating neutron 
stars. For simplicity, and in order to obtain a robust up- 
per limit, we obtain the gravitar limit under the following 
assumptions: all neutron stars are gravitars; all gravitars 
start emitting with the same initial frequency; this fre- 
quency is infinite; and the spectrum can be extended to 
arbitrarily low frequencies, as if signals had an infinite 
amount of time to evolve. Under these unrealistic as- 
sumptions, the energy spectrum can be approximated by 
Equation ( 77 ) , and the spectral function of the gravitar 



limit (using Equation (47)) becomes 



«gl(/) 



7T 2 I 

p c c 2 H 



h(z)£ 1 (z)dz. 



(103) 



Using the star formation rate of Equation ( 56 1 , the back- 
ground yielded by such an ensemble would produce a 
SNR (Equation p2"| ) of ~1.3 for aLIGO, after one year 
of observation. We use the gravitar limit as a reference 
in the following plots. 

We now justify that the background of rotating neu- 
tron stars cannot be larger than the gravitar limit. The 
spectral function in Equation ( 103 1 depends only on the 



rate n(z) and on the average moment of inertia /. The 
latter is well constrained by present neutron star equa- 
tions of state. The abundance and even the existence of 
gravitars is unknown, but certainly not all neutron stars 
are gravitars, so the rate of gravitars must certainly be 
smaller than h(z). Equation (103) is obtained by assum- 
ing that fc = 0. If the cut frequency were not zero, at 
frequencies lower than fc the spectral function would be 
proportional to / 4 , reaching its maximum around the cut 
frequency. There is hence no choice of the parameters B, 
e, and a, and there is no frequency at which the spec- 



tral function can be larger than ( 103 1, as long as the rate 
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FIG. 4. Spectral function, versus observed frequency, of the 
total and unresolvable backgrounds produced under the as- 
sumption that all rotating neutron stars are gravitars, i.e. 
neutron stars which spin down primarily by emitting gravi- 
tational waves. The gravitar limit is also shown, as a robust 
upper limit of the background from rotating neutron stars. 
The background is calculated by using three different initial 
frequency distributions (see Section IVE I. The unresolvable 
background is calculated with Ao = 1 and Af = lyr -1 . 



FIG. 5. Spectral function, versus observed frequency, of the 
total background produced by rotating neutron stars, assum- 
ing the magnetic field distribution of Equation ( 88 1 (which 
is taken from [77]) and the ellipticity distribution of Equa- 
tion (891 (from 15 ). Three initial frequency distributions 
are used (see Section |IVE[ ) . None of the present or planned 
gravitational wave detectors is sensitive enough to observe 
such a background. 



and the momentum of inertia (as well as the cosmological 
parameters) remain unchanged. 

If all rotating neutron stars were gravitars, the back- 
ground they would produce would be different than the 
gravitar limit. First, their initial frequency is finite, and 
second, they had a finite amount of time to evolve, so 
they cannot emit at arbitrarily low frequencies. In Fig- 
ure |4j besides the gravitar limit, we show the background 
that would be produced if all rotating neutron stars were 
gravitars, assuming the three initial frequency distribu- 
tions considered in Section IIV El The curves for Distri- 
bution are obtained by using Equation (43), whereas 



those for Distributions 1 and 2 are obtained by evaluat- 
ing Equation ( 33 ) . In all cases, the star formation rate 



is the one in Equation ( 56 1 , the magnetic field is approx- 
imately zero, and the ellipticity follows the distribution 
given in Equation ( 89 ) . The obtained background turns 



out to be almost entirely unresolvable. 

The SNR produced by the total background, assum- 
ing Distribution 0, is of 0.64, 6.6xl0 2 , and 3.5xl0 2 , for 
aLIGO, ETB, and ETD, respectively (assuming one year 
of observation). With Distribution 1, these numbers are 
0.56, 4.1 xlO 2 , and 1.8xl0 2 . Finally, with Distribution 
2, the values of SNR are 4.0xl0~ 3 , 2.0, and 0.81. We 
can thus claim that aLIGO is not sensitive enough to 
either detect the background of rotating neutron stars, 
or to set upper limits on the fraction of neutron stars 
that are gravitars. The SNR for BBO and DECIGO is 
in all cases much smaller than 1. Obviously, not all neu- 
tron stars are gravitars; if only a certain fraction of the 
population of neutron stars were gravitars, the values of 
the spectral function in Figure |4j as well as the values of 
SNR, would be multiplied by that fraction. Hence, if only 
1% of neutron stars were gravitars, they would produce 



a background that could be detected by ETB with SNR 
6.6, assuming Distribution 0, and 4.1, assuming Distri- 
bution 1. On the other hand, with Distribution 2, even 
if all neutron stars were gravitars the detection statistics 
of all detectors are below the detection threshold. 



B. A more realistic expectation 

In Figure [5] we show the spectral function of the total 
background of rotating neutron stars, calculated by as- 
suming the magnetic field distribution of Equation ( 88 1 , 



and the ellipticity distribution of (89). The star forma- 
tion rate used is the one of Equation ( 56 1 . The obtained 
SNR for such a background is much smaller than 1, for 
all detectors considered. 



C. Upper limit for magnetars 

There are two facts that make the detection of the 
magnetars background difficult: First, since they have a 
larger magnetic field, the electromagnetic emission dom- 
inates over the gravitational wave emission. Second, the 
fraction of magnetars is believed to be of order 10% of the 
total population of neutron stars [3H [79] . On the other 
hand, large magnetic fields can deform a neutron star 
|14) . increasing its ellipticity. If the ellipticity is large 
enough, the contribution of gravitational waves can be 
important. Furthermore, the larger the magnetic field, 
the faster systems evolve towards lower frequencies, en- 
tering the band of highly sensitive detectors like BBO 
and DECIGO. We now investigate if the background of 
magnetars has a good chance to be detected. 
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FIG. 6. Spectral function, versus observed frequency, of the 
total and unresolvable backgrounds produced by magnetars. 
The rate of magnetars is assumed to be 10% of the one of 
neutron stars, the average magnetic field is B = 10 10 T, and 
the average ellipticity, e = —6.4 x 10 -4 . This corresponds 
to the TD model described in [31], which is the most op- 
timistic model (regarding the detection) considered in that 
paper. Other models predict levels of background several or- 
ders of magnitude lower. The total background plotted here 
can thus be considered an upper limit for the background of 
magnetars. 

In [34], different models for the population of magne- 
tars are compared. One of the models, with a domi- 
nating toroidal magnetic field, produces a gravitational 
wave background that can be detected by ET. This model 
predicts a poloidal magnetic field of B = 10 10 T and an 
ellipticity e = —6.4 x 10~ 4 . In Figure [6] we show the 
background produced with this model, assuming that all 
magnetars start emitting gravitational waves of the same 
frequency, / max [82]. Other models in [34] predict levels 
of background orders of magnitude lower than the one 
in Figure [6] The total background shown in this figure 
can thus be considered an optimistic upper limit for the 
background of magnetars. 

The SNR with which the total background of Figure [6] 
would be detected by ETB and ETD is of 14 and 11, re- 
spectively (for one year of observation time). For aLIGO, 
as well as for BBO and DECIGO, the values of SNR are 
negligible. On the other hand, the SNR of the unre- 
solvable background is lower than 10~ 2 for all detectors. 
This means that the magnetars background is resolvable 
in the band of ET. 



D. Detection prospects 

In Section |"V A| we have claimed that the background 
produced by gravitars could be detected even if they con- 
stituted only a 1% of the neutron star population. Nev- 
ertheless, the existence of gravitars is questionable. We 
now show that similar detection claims can be achieved 
for certain (plausible) choices of B and e. 

In Figures [7] and [8] the total and unresolvable back- 
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FIG. 7. Spectral function, versus observed frequency, of the 
total background produced if all rotating neutron stars had 
the same magnetic field and ellipticity. Different line types 
correspond to different magnetic fields (as the legend of the 
upper plot describes). Each plot corresponds to a certain 
ellipticity. 
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FIG. 8. Same as Figure [7] but for the unresolvable part of 
the background. 



grounds are plotted, respectively, assuming that all ro- 
tating neutron stars have the same magnetic field and 
ellipticity. These plots show that the spectral function is 
larger for smaller magnetic fields and larger ellipticities, 
as it was expected. Furthermore, with larger magnetic 
fields, lower frequencies are achieved, and a bigger part 
of the background becomes resolvable. 

In Figures |9l [l0l[TH and[l2j we plot the SNR obtained 
by cross-correlating two interferometers of ETB, ETD, 
BBO or DECIGO, respectively, assuming one year of ob- 
servation. To obtain the corresponding SNR for an inter- 
val of observation time T Q b s , one can just multiply those 
values by -\/T bs/[l yr]. Each point on each curve corre- 
sponds to one realization of the background, in which all 
neutron stars have the same magnetic field and the same 
ellipticity. We know, of course, that not all neutron stars 
are equal. However, these plots are interesting for the 
following reason: the SNR (as well as the spectral func- 
tion) is proportional to the rate. So, all neutron stars 
may not have the same certain values of B and e, but if 
only a given fraction does, the SNR produced would be 
the one of the plots, multiplied by that fraction. This 
allows us to draw a few interesting conclusions. 

If Distribution 2 accurately describes the initial distri- 
bution of frequencies, the detection of the background of 
rotating neutron stars seems unlikely; instead, if Distri- 
butions or 1 are accurate, the chances of detection (by 
ET) are high. For example, we see in Figure [9] (assuming 
Distribution 1) that, if all neutron stars have an elliptic- 
ity of e = 10 -6 and a magnetic field of B = 10 7 T, the 
obtained SNR is 1.0 x 10 2 for ETB; then, if at least a 
few percent of neutron stars have ellipticities larger than 
10~ 6 and magnetic fields smaller than 10 7 T, ETB would 
detect the produced background with SNR of at least a 
few. Suppose now that neutron stars cannot have ellip- 
ticities larger than 10~ 7 . Even in this case, SNR of a few 
would be obtained for ETB if only a few percent of the 
rotating neutron stars have magnetic fields lower than 
10 6 T. 



As Figures 11 and 12 show, the SNR calculated for 
BBO and DECIGO reach relevant values for magnetic 
fields smaller than ~ 10 7 T, and, in Figure [7J we see 
that the spectral function, for such a magnetic field, 
has support only at frequencies larger than «lHz. The 
main contribution to the SNR of BBO and DECIGO 
thus comes from frequencies between 1Hz and 10 Hz. In 
[83] . the overlap reduction function is calculated for dif- 
ferent configurations of the spacecraft constellations of 
BBO and DECIGO. Almost all configurations produce 
an overlap reduction function close to zero between 1 Hz 
and 10 Hz. The assumption made in Section [II C| of an 
overlap reduction function equal to one is therefore very 
crude. The SNR obtained with a more realistic overlap 
reduction function would reasonably be much lower. 



In Figures 13 and [14| we show the same SNR calcu- 
lations as in Figures |9j and [Toj but for the unresolvable 
part of the background. By comparison, one can con- 
clude that, for ellipticities smaller than e ~ 10~ 4 and 
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FIG. 9. Signal-to-noise ratio produced by the total back- 
ground of rotating neutron stars, assuming that all of them 
have the same magnetic field and ellipticity. Each curve cor- 
responds to a certain magnetic field (specified on top of each 
curve), and each point on a curve corresponds to a certain 
ellipticity (specified on the horizontal axis). These values of 
SNR are obtained by cross-correlating 1 year of data of two 
interferometers of ETB. Upper, middle, and lower plots are 
obtained using the initial frequency distribution called Distri- 
bution 0, 1, and 2, respectively, in Section [IV E| 



magnetic fields smaller than B 
is almost entirely unresolvable. 



10 8 T, the background 



FIG. 10. Same as Figure [9] calculated in this case for ETD. 



VI. DISCUSSION 

A. Comparison with previous work 

The background produced by magnetars has recently 
been calculated in |34) . assuming different models. Our 
upper limit is obtained by assuming the magnetic field 
and ellipticity of one of the models given in that paper: 
the one that predicts the largest spectral function. All 
other models in that work, as well as the models used in 
previous papers [33] [33] produce smaller levels of back- 
ground. 

In [S3], the gravitational wave background is calcu- 
lated, assuming a population of magnetars that could fit 
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FIG. 11. Same as Figure [9] but calculated for two interfer- 
ometers of BBO. The observation time assumed is also of 1 
year. 



FIG. 12. Same as Figure [9] but for two interferometers of 
DECIGO. The observation time assumed is also of 1 year. 



the ultrahigh energy cosmic ray spectrum. The most op- 
timistic of the expectations for the spectral function in 
that paper is, in a certain range of frequencies, a factor of 
ss 8 larger than the upper limit calculated by us in Sec- 
tion |V C| That expectation of the background is said to 
possibly reach the sensitivities of BBO and DECIGO, but 
not that of ET; however, no calculation of SNR is per- 
formed in that paper. We now assert that, performing 
the SNR calculations, the claim is the opposite. The SNR 
that our magnetar upper limit would produce on BBO 
and DECIGO (assuming an overlap reduction function 
equal to one, which is already too optimistic), after one 
year of observation time, is of the order of 10~ 3 ; these 



values are too low to claim a possible detection, even if 
multiplied by that factor of ~ 8. However, the upper 
limit of [53] would produce an SNR of as 5 and a 2 on 
ETB and ETD, respectively. Therefore, one can conclude 
that the upper limit of magnetars (either with the esti- 
mate of [83] or with ours) is out of the reach of BBO and 
DECIGO, but could be detected by ET. 

All the papers mentioned in this section use the so- 
called duty cycle to account for the statistical properties 
of the background. As commented in Section V.F.3 of 
|18j . the overlap function (which is a generalization of 
the duty cycle) is the right tool to quantify the resolv- 
ability of the background. Moreover, the duty cycle can 
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FIG. 13. Same as Figure [9] but for the unresolvable part 
of the background (using Ao = 1 and A/ = lyr^ 1 ). These 
values of SNR are obtained by cross-correlating 1 year of data 
of two interferometers of ETB. 

only be used for short (burst-like) signals, not for long 
signals, like the ones produced by rotating neutron stars. 
In Figures [6j [l3j and [Til one sees that, even if having a 
duty cycle much larger than 1 (as reported in [M]), the 
background produced by magnetars is resolvable in the 
band of ET. 



B. On the formation rate 

The spectral function turns out to be rather insensitive 
to the shape of the rate. In this section we compare 



FIG. 14. Same as Figure [Tj) but for ETD. 

the spectral function obtained using the different star 
formation rates piz) of and a star formation rate 

that has the same value over all cosmological epochs. 

Let us first obtain a reasonable value for the constant 
rate n(z) — R. Given one star formation rate f>i(z), one 
can calculate its average value over all redshifts, 

f " 1 ^ p i (z)dz 

(Pi) = J " rJl , ■ (104) 

Jo dz 

Considering the N = 11 star formation rates from fS7\ - 
[74] . the mean value of (pi) is 

i N 

W = N £te> = °- 1OM0 ^ MPC " 3 ■ (1 ° 5) 
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FIG. 15. Spectral function, versus observed frequency, of 
the total background produced by the ensemble of rotating 
neutron stars, using the three initial freque ncy distributions 
(Distribution 0, 1, and 2, of Section IVEl. All systems are 
assumed to have B 



l(fT, and t = 10"'. The light-gray 
curves are obtained with different star formation rates [67H74] , 
whereas the black dashed lines use the constant star forma- 
tion rate of Equation (106 1. The black solid line is obtained 



with the approximate formula of Equation ( A7 I 



Finally, replacing this star formation rate in Equation 
( 57 1 , one obtains the value of the constant rate R, 



R = X(p) = 5x 10~ 4 yr" 1 Mpc 



-3 



(106) 



This is the approximate number of neutron stars formed 
per unit emitted interval of time per unit comoving vol- 
ume, at any cosmological epoch. 

In Figure |15| we show the spectral function, with the 



three initial frequency distributions of Section IV E ob- 
tained with constant and non-constant rates. One sees 
that the spectral function using different rates differs just 
by a constant factor at almost all frequencies; only at fre- 
quencies close to the maximum one, the spectral function 
presents different shapes depending on the specific rate 
assumed. Given the lack of information on the popu- 
lation statistics (which is evident when comparing the 
results of different initial frequency distributions), the 
choice of one or another rate is irrelevant. A similar con- 
clusion is drawn for binary systems in Section V.A of 



that evolves over long time scales. We also present an ex- 
pression for the overlap function Af(f, A/, z), that gives 
the number of signals with observed frequency / per fre- 
quency bin A/, with redshifts smaller than z. This over- 
lap function (in Equation (16)) is more general than the 
one introduced in [15] . 

A more detailed expression of the spectral function is 
given in Equation (33), where the dependence on the ini- 



tial frequency distribution of the ensemble is explicitly 
shown. Similarly, a more explicit formula for the over- 
lap function is presented in Equation (34) (or, in a more 



convenient form, in Equation (35)). Assuming that all 
systems start emitting at the same frequency, the over- 
lap function and the spectral function become the ones 
of Equations (42) (or, more conveniently, ([44]) ) , and (43), 
respectively. Equation ( 46 1 gives the overlap function 



calculated with the additional assumption that the evo- 
lution of the systems is short compared to cosmological 
time scales. This formula coincides with the definition 
given in [T5J. The spectral function is then obtained un- 
der the same assumptions (Equation (47)). If one applies 
the latter formula to calculate the total background (with 
Afo = 0), one obtains Equation (48), which is the expres- 
sion that usually appears in the literature. 

In Section [TV] we have described a simple but useful 
model for the energy and frequency evolution of rotating 
neutron stars. Figures [2] and [3] show the collection of pos- 
sible redshifts and observed frequencies of gravitars and 
of normal pulsars, respectively, assuming in both cases 
that all systems have the same ellipticity, e = 10~ 7 . In 
these plots one can see the frequency range in which the 
usual assumption of short-lived systems (which has not 
been adopted for our results) is inaccurate. 

In Appendices [A] and [E] we propose simple approxi- 
mate formulas for the spectral function of rotating neu- 
tron stars and, for completion, for binary systems too. 
Given the uncertainties in some astrophysical parame- 
ters, like the rate and the initial frequency distribution, 
these approximate formulas can be used as a reason- 
ably good estimation for the levels of contemporary back- 
grounds. 



B. Regarding the detection of the background of 
rotating neutron stars 



VII. SUMMARY AND CONCLUSIONS 



A. Regarding the calculation of the spectral 
function 



We have derived a formula for the spectral function 
fi(/, Af,Afo), i.e., the density per logarithmic frequency 
interval (in units of critical density), around the observed 
frequency /, of a background made of signals that pro- 
duce an overlap of Afo signals per frequency bin A/. This 
formula (in Equation (18)) can be used for a population 



The three initial spin period distributions considered 



in the calculations (see Section IV E), lead to very dis- 
parate results. The first one (called Distribution 0), as- 
sumes that all systems are formed with the same initial 
frequency. The other distributions (called Distributions 
1 and 2) are taken from the literature [35] [77] . If Distri- 
bution 2 is accurate, the detection of the background by 
present and planned detectors can be discarded. 

In Section [V A| we have justified that there is a robust 
upper limit (the gravitar limit, in Equation (103), or, 



more simply in ( A5 )) on the level of background produced 
by rotating neutron stars. In Appendix O we point out 
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an interesting characteristic of the background obtained 
under the unrealistic assumptions of the gravitar limit: 
its characteristic amplitude h c is independent of the ellip- 
ticity of the systems and of their spin frequency. A simi- 
lar feature was predicted by Blandford, for the expected 
gravitational wave amplitude of the nearest system of a 
uniform galactic population of gravitars. This feature 
disappears with more realistic models, as it occurs with 
Blandford's argument. 

In Figure |4j the spectral functions of the total and the 
unresolvable background are calculated, under the un- 
realistic assumption that all neutron stars are gravitars. 
One sees that the background is almost entirely unre- 
solvable. Such a background, assuming Distributions 
and 1, can be detected with ETB and ETD, but not with 
aLIGO, by using the typical cross-correlation method. If 
at least 1% of neutron stars behave like gravitars, they 
will produce an unresolvable background that can be de- 
tected by ET. 

Using a magnetic field and an ellipticity distribution 
(from [T21 177]), a reasonable level of background of ro- 
tating neutron stars is obtained (see Figure [5]); it is be- 
low the detection capabilities of any existing or planned 
instrument. 

We have obtained an upper limit on the total back- 
ground of magnetars (Figure [6]) , using one of the models 
presented in [33]. This background can only be detected 
by ET. However, other models predict levels of back- 
ground several orders of magnitude lower. Hence, we 
claim that magnetars are not the most promising rotating 
neutron stars, regarding the detection of the background. 

With Figures [7] and [8] one can get an idea of how the 
amplitude and the resolvability of the background depend 
on the values of the magnetic field and ellipticity (assum- 
ing that these values are equal for all neutron stars). 

Figures |9j [lOj [IT] and [12] summarize the prospects of 
ETB, ETD, BBO, and DECIGO, respectively, to detect 
the total background of rotating neutron stars. The val- 
ues of SNR in these plots are calculated for the cross- 
correlation of the data of two interferometers during 1 
year (to obtain the values after 3 or 5 years, it is enough 
to multiply the vertical axis by or y5) ■ These graphs 
are useful because the SNR is proportional to the fraction 
of stars that are neutron stars. Suppose that some given 
values of magnetic field B and ellipticity e are associated 
with SNR equal to S in the plots; then, if only a fraction 
x of all neutron stars have B and e, they will produce a 
background observed with SNR equal to XX S. For exam- 
ple, one can conclude from Figure [9] that, if at least a few 
percent of neutron stars have e > 10 -6 and B < 10 7 T, 
the background would be observed by ETB with SNR of 
a few. Another conclusion from Figure [9] is that, if the 
maximum ellipticity of neutron stars is of e = 10 , ETB 
will still observe a background of SNR of a few, if just a 
few percent of neutron stars have B < 10 6 T. 

We point out that the SNR values of BBO and DE- 
CIGO are obtained by assuming an overlap reduction 
function (see Section II C I identically equal to 1. This is 



quite inaccurate between 1Hz and 10 Hz (see [H3], where 
the overlap reduction function is calculated for different 
configurations of the spacecrafts) , which is the frequency 
interval where the background contributes the most to 
the SNR of BBO and DECIGO. The detection prospects 
of BBO and DECIGO should therefore be smaller than 
what Figures [IT] and [T2] suggest. 

Figures [13] and [TJJ are analogous to [9] and [TU] respec- 
tively, but for the unresolvable part of the background. 
They are included to show that the total and the unre- 
solvable backgrounds are identical for all configurations 
with e < 10~ 4 and B < 10 8 T; on the other hand, the 
background produced by magnetars (with larger mag- 
netic fields) is mostly resolvable. 

This work, together with [TS] , covers two of the most 
promising sources of contemporary gravitational wave 
background. If the most reasonable estimate of the back- 
ground (in Section VB) is accurate, or if rotating neu- 



tron stars form with initial spin frequencies well described 
by Distribution 2 (in Equation (85)), then we can con- 
clude that ground-based detectors operate in a frequency 
window that is free of contemporary unresolvable back- 
ground from binaries and rotating neutron stars. How- 
ever, if at least a few percent of neutron stars behave 
as gravitars, or if any of the configurations proposed in 
Section |VD| that produce high SNR values is in good 
agreement with the real neutron star population, an un- 
resolvable background of rotating neutron stars can ob- 
scure the searches of ET for other sources of background. 
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Appendix A: Simple formulas for the background of 
rotating neutron stars 

In this section we present some approximate formulas 
for the spectral function of the total and the unresolvable 
background of rotating neutron stars. They are useful to 
obtain a simple estimate of the background; neverthe- 
less, these formulas were not used to obtain the results 
of Section [V] 

We first need to define the constants 



C gr — 



S~ 1 (z)dz 1.9, 



Cmd — 



[1 + z] 2 £- 1 (z)dz = 



2i£ ( z max ) 



(Al) 

19.0, 
(A2) 
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Cmd — 

and 

C gr = 



Jo Uo 



£' 1 {z)dz « 0.23, 
(A3) 

£- 1 (z)dz ps 0.042, 
(A4) 

which depend only on cosmological parameters. 

A formula for the spectral function of the gravitar limit 
can be very easily obtained. Using a constant rate h(z) — 
R, Equation ( 103 ) becomes 



[1+*] 



£~ x {z!)dz l 



A value for i? can be found in Section IVIBI 

In the following, we assume a fixed initial frequency 
/max, and obtain the spectral function by solving Equa- 
tion (47) for rotating neutron stars. Moreover, we as- 
sume a constant rate n(z) = R, and perform the approx- 



imation introduced in Section IV C of distinguishing md- 
and gr-ranges. The spectral function of the background 
of rotating neutron stars, under these assumptions, has 
an analytical form. To account for the gr- and the md- 
regimes, we define the function 



«(/) 



^1R C gr / 2 



p c H c 2 



C md / 4 /< 



fc<f 
5fc 



(A6) 



where the cut frequency fc is the one in Equation ( 72 ) 
The total background can then be estimated using 



^total(Z) 




(A7) 



The maximum frequency is the one defined in Equation 
( 54 ) , and the minimum frequency is 



/min 



AOgr) fc<h(z&) 
/2 (fe+^) / 1 (z gr )</ c </ max 

fa(Zmd) /max < fc 



(A8) 



where fi{z), h{z), f 3 (z), z gr , and z md are given in 



Equations J94|, {95}, (|96j), ( |10l[ ), and ([102]), respectively. 
Equation JAM is an approximation to the minimum value 
of fiow(z), given by Equation (93); / m ; n is, therefore, de- 
fined as an observed frequency (unlike /max, that is an 
emitted frequency). 

The spectral function of the unresolvable background 
can be approximated by 



unresolvable 



(/) 



o 

^total(/) /] 





/</p .nun 

p,min — f — /p,max 
/p,max ^ f 



(A9) 



In this equation we have introduced the limiting frequen- 
cies (see Section III.E.4 of [H]); / p , m i n (/ P ,max) is the 
minimum (maximum) frequency above (below) which the 
unresolvable background is present. The limiting fre- 
quencies can be obtained from 



/p 



and 



min (x, f n 



mm x 



3/5 



2/5 



(A5) 

where 



(A10) 



X<fc 

i /maxj fc *^ X 

(All) 



8nc 6 RAf8 md C 



7^ lV3 



(A12) 



The upper case in (All) occurs when the unresolvable 



background is restricted solely to the md-range, and the 
lower case occurs when the unresolvable background is ei- 
ther restricted to the gr-range or partially in both ranges. 
Instead of / p ,min and / p , m ax, one can use / d , m in and / d ,max 
in Equation (A9); /d, m in (/d.max) is the minimum (max- 



imum) frequency above (below) which the unresolvable 
background dominates over the resolvable. The unresolv- 
able background is said to dominate over the resolvable 
when the spectral function of the former is larger than 
the spectral function of the latter. One can prove that 
/d.min ~ / P ,min, whereas /d.max is constrained to the in- 
terval 



*7p,max < /d,max < /, 



p,max • 



(A13) 



The factor F depends only on cosmological parameters, 
and is F ps 0.9 (both in the md and in the gr ranges). 
With Equation ( A13 ) we can conclude that, as soon as an 



unresolvable background appears at a certain frequency 
/ P ,max, that background dominates over the resolvable 
background below ss 0.9/ p max . 



Appendix B: Simple formulas for the background of 
binary systems 

To have a complete estimate of some of the most 
promising sources of contemporary backgrounds, we also 
give some simple approximate formulas regarding the 
background of stellar binary systems (those systems com- 
posed of neutron stars, white dwarfs, or stellar-mass 
black holes). These formulas are based on the calcula- 
tions of [IB]. 

The spectral function of the total background can be 



calculated again using Equation ( A7 1 , but with different 
definitions of ri(/), /max, and / m i n . For binary systems, 
we have 



_ #bm[G7r] 2 / 3 miro 2 Cbin ^2/3 



(Bl) 
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Here, mi and m-i are the masses of the two components of 
the binary, i?bin is the binary rate (summarized in Table 
I of [IB]), and 



Cbin — 



' '[1 + z]-^ 3 £- 1 {z)dz^ 0.75. 



(B2) 



The maximum frequency can be the frequency of the last 
stable orbit, which is given by 



where 

Cbin 



-8/3 



£- l {z')dz' 



£- 1 {z)dz « 0.12. 

(B7) 



In Equation (A9|, instead of / p , m m and / p , m ax, one can 
use / d ,min (which is approximately / p , min ) and / d ,ma X ; 



the latter also fulfills Equation (A13), but, in the case of 
binaries, the value of the F-factor is approximately 0.6. 



/max — 



mi + 7712 



(B3) 



For systems containing a white dwarf, a better choice of 
the maximum frequency is 



fWD 
' max 



'G[m\ + m%\ 
?: 2 [ri + r 2 ] 3 



(B4) 



where r\ and r 2 are the radii of the components; this 
frequency corresponds to a separation of the two compo- 
nents equal the sum of their radii. On the other hand, 
the minimum frequency is approximately given by 



256T(z max ) 7 r 8 / 3 G 5 / 3 m 1 m 2 



-i -3/8 



5c 5 [mi + m^] 1 / 3 



[1 ^max] 



(B5) 

where the function T{z) is the same function that has 
been used in the previous sections, defined in Equation 
( 25 1 . The minimum frequency for binary systems has 



been defined as an observed frequency (unlike in |18j ) in 
analogy to the minimum frequency defined for rotating 
neutron stars, in Equation (A8). 



To obtain the spectral function of the unresolvable 



background, one can once more use Equa tion (A9), with 



the definition of £l(f) given in Equation (Bl|, and with 
the limiting frequencies giv en b elow. First, / p , m i n ~ /minj 
which is given in Equation ( B5 ) . Second, the limiting fre- 
quency / P: max is approximately 



/ P ,max 



5A/C bin c 8 i? bin [m 1 +m 2 ] 1 / 3 
247r 5 / 3 G 5 / 3 m 1 m 2j ff 3 A/'o 



3/11 



) fn 



(B6) 



Appendix C: On the connexion between the gravitar 
limit and Blandford's argument 

Let us consider again the gravitar limit, described by 
the spectral function in Equation ( A5 1 . If a stochastic 



background was characterized by such a spectral func- 
tion, the characteristic amplitude of the gravitational 
wave spectrum |45) would be 



where 



(Cl) 



(C2) 



Using Equation (|A5j), we obtain a characteristic ampli- 
tude of the form 



h c 



2irGIRC s 



H c 2 



(C3) 



which does not depend either on the frequency or on the 
ellipticity of the systems. This fact reminds us Bland- 
ford's argument (first cited in [6], revisited in [85] and 
[16]). According to this argument, the expected gravita- 
tional wave amplitude of the nearest system of a uniform 
galactic population of gravitars, is independent of the el- 
lipticity and the spin frequency of the systems. Some of 



the assumptions required to obtain Equation ( C3 1 are dif- 
ferent from those of Blandford's argument. However, the 
result is similar: the expected gravitational wave charac- 
teristic amplitude of the background produced by a popu- 
lation of gravitars, uniformly distributed in the universe, 
is independent of the ellipticity and the spin frequency 
of the systems. 

The assumptions needed for Blandford's argument to 
hold, as well as those needed for h c not to be a func- 
tion of e and /, are quite unrealistic. Once one adopts 
more realistic assumptions on the galactic population (for 
example, that gravitars are not distributed on a two- 
dimensional disk), Blandford's argument vanishes [16] . 
Analogously, under more realistic assumptions on the en- 
semble of gravitars in the universe (for example, that they 
have a finite time to evolve, and a finite initial spin fre- 
quency), the characteristic amplitude in Equation (C3) 
depends on the ellipticity and on the frequency. 



24 



[1] B. Allen and J. D. Romano, Phys. Rev. D 59, 102001 
(1999). 

[2] M. Maggiore, Phys. Reports 331, 283 (2000). 
[3] J. M. Lattimer and M. Prakash, Science 304, 536 (2004). 
[4] A. J. Deutsch, Annales d'Astrophysique 18, 1 (1955). 
[5] F. Pacini, Nature 219, 145 (1968). 

[6] Section 9.4.2 (b) of S. W. Hawking and W. Israel, 300 
Years of Gravitation (Cambridge University Press, Cam- 
bridge, UK, 1987). 

[7] R. Prix, in Neutron Stars and Pulsars, Astrophysics 
and Space Science Library, Vol. 357, edited by 
W. Becker (Springer Berlin Heidelberg, 2009) pp. 651- 
685, 10.1007/978-3-540-76965-l_24. 

[8] D. Lorimer and M. Kramer, Handbook of Pulsar As- 
tronomy (Cambridge University Press, Cambridge, UK, 
2005). 

[9] D. Lorimer, Living Rev. Relativity 11, 8 (2008). 
[10] C. Thompson and R. C. Duncan, Astrophys. J. 408, 194 
(1993). 

[11] C. Thompson and R. C. Duncan, Mon. Not. R. Astron. 
Soc. 275, 255 (1995). 

[12] C. Thompson, M. Lyutikov, and S. R. Kulkarni, Astro- 
phys. J 574, 332 (2002). 

[13] A. K. Harding and D. Lai, Rept.Prog.Phys. 69, 2631 

(2006) . 

[14] C. Cutler, Phys. Rev. D 66, 084025 (2002). 
[15] C. Palomba, Mon. Not. R. Astron. Soc. 359, 1150 (2005). 
[16] B. Knispel and B. Allen, Phys. Rev. D 78, 044031 (2008). 
[17] S. B. Popov et al, Astrophys. J 530, 896 (2000). 
[18] P. A. Rosado, Phys. Rev. D 84, 084004 (2011). 
[19] B. J. Owen et al., Phys. Rev. D 58, 084020 (1998). 
[20] X.-J. Zhu, X.-L. Fan, and Z.-H. Zhu, Astrophys. J 729, 
59 (2011). 

[21] L. Barack and C. Cutler, Phys. Rev. D 70, 122002 (2004). 
[22] P. Amaro-Seoane et al, Clas. Quantum Grav. 24, R113 

(2007) . 

[23] A. Buonanno et al, Phys. Rev. D 72, 084001 (2005). 

[24] S. Marassi, R. Schneider, and V. Ferrari, Mon. Not. R. 
Astron. Soc. 398, 293 (2009). 

[25] I. Kowalska, T. Bulik, and K. Belczynski, As- 
tron. Astrophys. 541, A120 (2012). 

[26] B. Allen, arXiv: gr-qc/9604033v3 (1996). 

[27] A. Buonanno, arXiv: gr-qc/0303085v2 (2004). 

[28] C. Cutler and J. Harms, Phys. Rev. D 73, 042001 (2006). 

[29] C. Cutler and D. E. Holz, Phys. Rev. D 80, 104009 

(2009) . 

[30] K. Yagi and N. Seto, Phys. Rev. D 83, 044011 (2011). 
[31] T. Regimbau and J. A. de Freitas Pacheco, Astron. As- 
trophys. 376, 381 (2001). 
[32] T. Regimbau and J. A. de Freitas Pacheco, Astrophys. J 

642, 455 (2006). 

[33] T. Regimbau and V. Mandic, Clas. Quantum Grav. 25, 
184018 (2008). 

[34] S. Marassi et al, Mon. Not. R. Astron. Soc. 411, 2549 
(2011). 

[35] C.-A. Faucher-Giguere and V. M. Kaspi, Astrophys. J. 

643, 332 (2006). 

[36] M. Punturo et al, Clas. Quantum Grav. 27, 194002 

(2010) . 

[37] H. Grote, Clas. Quantum Grav. 27, 084003 (2010). 
[38] T. Accadia and B. L. Swinkels, Clas. Quantum Grav. 27, 



084002 (2010). 

[39] LIGO Scientific Collaboration, RPP 72, 076901 (2009). 
[40] G. M. Harry, Clas. Quantum Grav. 27, 084006 (2010). 
[41] S. Kawamura et al, Clas. Quantum Grav. 28, 094011 
(2011). 

[42] A. G. Riess et al, Astrophys. J 699, 539 (2009). 

[43] A. G. Riess et al., Astrophys. J 730, 119 (2011). 

[44] N. Jarosik et al, Astrophys. J. Suppl. Ser. 192, 14 (2011). 

[45] E. S. Phinney, arXiv: astro-ph/0108028vl (2001). 

[46] A more thorough definition of the resolvability, that takes 

into account the difference in amplitude of the signals, 

can be the subject of a future work. 
[47] T. Regimbau et al., arXiv: gr-qc/1201.3563vl (2012). 
[48] M. Maggiore, Gravitational Waves Volume 1: Theory 

and Experiments (Oxford University Press, New York, 

USA, 2008). 

[49] B. F. Schutz, pas. Quantum Grav. 28, 125023 (2011)"] 

[50] L. S. Finn, S. L. Larson, and J. D. Romano, Phys. Rev. 
D 79, 062003 (2009). 

[51] Taken from LIGO Scientific Collaboration, Advanced 
LIGO anticipated sensitivity curves, Tech. Rep. LIGO- 
T1000414-vl3 (LIGO Scientific Collaboration, 2010). 

[52] The spectral strain sensitivities of ETB and ETD were 
kindly provided by Tania Regimbau in a private commu- 
nication. 

[53] The spectral strain sensitivities of DECIGO and BBO 
are obtained by using the fitting formulas given in 
A. Nishizawa et al, Phys. Rev. D 85, 044047 (2012). 

[54] P. J. E. Peebles, Principles of Physical Cosmology, 
Princeton series in physics (Princeton University Press, 
Princeton, USA, 1993). 

[55] By imposing £ = 0, the function T(z) gives the age of 
the universe at the instant when the waves of redshift z 
were emitted, as in Equation (13.20) of |54| . 

[56] B. Knispel, Pulsar Discoveries by Volunteer Distributed 
Computing, Ph.D. thesis, Leibniz Universitat Hannover 
(2011). 

[57] One can prove it by partially differentiating Equation 

( |81[ ) with respect to t e . 
[58] T. Regimbau, Research in Astron. Astrophys. 11, 
369 (2011). See corrected version in arXiv: astro- 
ph/1101.2762v3. 



[59] In this regard, one should read Section IV G the plots in 
this section, for instance the ones in Figure d3j, can be 
qualitatively compared with Figures (2) and (4) of |18| . 

[60] J. M. Lattimer, New Astronomy Reviews 54, 101 (2010). 

[61] J. E. Trumper, Progress in Particle and Nuclear Physics 
66, 674 (2011). 

[62] J. D. Brown, Phys.Rev.D 62, 084024 (2000) 



[63] See Section 3.5 of D. R. Lorimer, High-Energy Emis- 
| sion fro m Pulsars and their Systems^ edited by D. Torres 

and N. Rea, Astrophysics and Space Science Proceedings 

(Springer Verlag, Berlin Heidelberg, 2011) pp. 21-36, and 

references therein. 
[64] ATNF catalogue: R. N. Manchester et al, Astron.J. 129, 

1993 (2005). 

[65] D. C. Backer et al, Nature 300, 615 (1982). 

[66] See Section 7.4 of S. L. Shapiro and S. A. Teukolsky, 
Black Holes, White Dwarfs, and Neutron Stars: The 
Physics of Compact Objects, Physics textbook (Wiley, 
New York, USA, 1983). 



25 



[67] P. Madau, M. D. Valle, and N. Panagia, Mon. Not. R. 

Astron. Soc. 297, L17 (1998). 
[68] C. Porciani and P. Madau, Astrophys. J 548, 522 (2001). 
[69] 2dFGRS Team, Mon. Not. R. Astron. Soc. 326, 255 

(2001). 

[70] L. G. Strolger et al, Astrophys. J 613, 200 (2004). 

[71] A. M. Hopkins and J. F. Beacom, Astrophys. J 651, 142 

(2006) . 

[72] K. Nagamine et al, Astrophys. J 653, 881 (2006). 

[73] M. A. Fardal et al, Mon. Not. R. Astron. Soc. 379, 985 

(2007) . 

[74] S. M. Wilkins, N. Trentham, and A. M. Hopkins, Mon. 

Not. R. Astron. Soc. 385, 687 (2008). 
[75] E. E. Salpeter, Astrophys. J 121, 161 (1955). 
[76] S. Johnston and D. Galloway, Mon. Not. R. Astron. Soc. 

306, L50 (1999). 
[77] Z. Arzoumanian, D. F. Chernoff, and J. M. Cordes, As- 



trophys. J. 568, 289 (2002). 
[78] T. Regimbau an d J. A. d e Freitas Pacheco, Astron. As- 



[79] 



trophys. 359, 242 (2000) 
B 



Mon. Not. R. Astron. Soc. 401, 2675 



S. B. Popov et al. 
(2010). 

[80] R. Perna et al, Mon. Not. R. Astron. Soc. 384, 1638 
(2008). 

[81] P. L. Gonthier et al, AIP Conf. Proc. 1357, 245 (2011). 

[82] Some studies predict that magnetars are formed with fast 
initial spins [TO]. Since we are interested in obtaining an 
upper limit, we assume that all magnetars start emitting 

^1 /max - 

[83] A. Nishizawa, A. Taruya, and S. Kawamura, Phys. Rev. 

D 81, 104043 (2010). 
[84] K. Kotera, Phys. Rev. D 84, 023002 (2011). 
[85] LIGO Scientific Collaboration, Phys. Rev. D 76, 082001 

(2007). 



